Homodyne detection of matter-wave fields 
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A scheme is discussed that allows one for performing homodyne detection of the matter-wave 
field of ultracold bosonic atoms. It is based on a pump-probe lasers setup, that both illuminates 
a Bose-Einstein condensate, acting as reference system, and a second ultracold gas, composed by 
the same atoms but in a quantum phase to determine. Photon scattering outcouples atoms from 
both systems, which then propagate freely. Under appropriate conditions, when the same photon 
can either be scattered by the Bose-Einstein condensate or by the other quantum gas, both flux 
of outcoupled atoms and scattered photons exhibit oscillations, whose amplitude is proportional to 
the condensate fraction of the quantum gas. The setup can be extended to measure the first-order 
correlation function of a quantum gas. The dynamics here discussed make use of the entanglement 
between atoms and photons which is established by the scattering process in order to access detailed 
information on the quantum state of matter. 



I. INTRODUCTION 

The measurement technique which is mostly employed 
in experiments with ultracold atoms is based on time of 
flight [J Q • Used in different setups it allows one to deter- 
mine the quantum state of matter ■ Further exciting 
perspectives have been opened by the recent demonstra- 
tion of in-situ imaging of atoms in optical lattices 
In addition, the remarkable progress in coupling ultracold 
atomic gases with high-finesse optical resonators [lol - fl3j 
has generated a renewed interest in revealing the corre- 
lation functions of matter by photodctcction of the scat- 
tered light. This progress opens, amongst several others, 
the possibility to monitor properties of the quantum state 
of matter in a non-destructive way fl4l - |l~7j . 

Photons emitted by Rayleigh scattering, and in gen- 
eral in a pump-probe type of experiment, have been used 
to characterize the quantum state of atomic gases [f8l - 
[22j . They deliver information on the structure form fac- 
tor of the quantum gas and thus on the density and 
density-density correlations 0, [23 - [2^ | . Proposals for op- 
tical detection of certain quantum states of matter have 
appeared in the literature [231 - 4271 ] - 

In this article we discuss a setup in which detection of 
the photons scattered by a quantum gas permits one to 
determine the mean value of the atomic field. This setup 
constitutes a matter-wave analogon of homodyne detec- 
tion of light fields [28| , and specifically allows one to de- 
termine the condensate fraction of a quantum gas by 
photodetection. Furthermore we show that this scheme 
can be extended to determine the first-order correlation 
function of an atomic gas when the two illuminated re- 
gions belong to the same quantum gas. 

The setup for determining the condensate fraction 
is sketched in Fig. Q] and is based on an interferom- 
eter for Bose-Einstein Condensates (BEC) realized in 
Refs. [13, HJ. We establish a direct connection be- 
tween our theoretical model and the experimental setup 
in Ref. 130, Sll and then focus on a scheme in which the 



second system is not necessarily a BEC. We then show 
that this setup allows one to determine the condensate 
fraction of the second system by measuring the flux of the 
scattered photon. The dynamics is based on first creating 
entanglement between the scattered photon and the scat- 
tering atom of the quantum gas and then on performing 
an operation similar to a quantum eraser [32| , thereby al- 
lowing one to measure the correlation functions of matter 
in the flux of the scattered light. 

This article is organized as follows. In Sec. |TI]the setup 
and the theoretical model are introduced. In Sec. IIIII 
a formal connection is derived between the photon flux 
and the correlation functions of the scattering atoms. In 
Sec. IIVI the photon flux is evaluated when the scatter- 
ing systems are two Bose-Einstein condensates at zero 
temperature and at different temperatures. In Sec. [V] it 
is discussed how a similar setup can be used to measure 
the first-order correlation function. The conclusions are 
drawn in Sec. IVII The appendices report details of the 
calculations presented in Sees. [Til IIIII and IIVI 



II. THE MODEL 



A realization of the setup we consider in this paper is 
depicted in Fig. [T] Here, ultracold bosonic atoms of mass 
m are initially prepared in the stable electronic state |1) 
and confined by the (state-dependent) potential Vf(r). 
The atoms that are illuminated by the lasers are coher- 
ently coupled to the electronic state |2), which is stable 
and not confined: the outcoupled atoms propagate freely 
in space. 

In this section we introduce the Hamiltonian and the 
physical quantities which are at the basis of our analysis. 
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FIG. 1: Setup for performing homodyne detection of matter- 
wave fields. A gas of identical bosonic atoms is confined in 
two distinguishable regions of space. In the left region it forms 
a Bose-Einstein condensate that acts as reference system for 
measuring the condensate fraction of the right system. A 
pump-probe laser setup outcouples the atoms from both wells. 
The photon flux of the probe beam is constituted by the pho- 
ton scattered in the outcoupling process and exhibits oscilla- 
tions whose amplitude depend on the condensate fraction of 
the right system, and which vanish when this is zero. The 
first-order correlation function of a quantum gas can be de- 
termined in an extension of this setup and is discussed in Sec. 
FVl Details on the parameters are given in the text. 



A. Hamiltonian 

We denote by H the Hamiltonian governing the dy- 
namics of photons and atoms, which we decompose into 
the sum 



as Hq — Hi + H2 + H±2, where 

, r 2 ^ 2 

Hi 4— 



b=i,2} = J drVj(r) + Vj(r) + hu a ^ ^(r) 

(2) 



ii 
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drip] (r)t/>j (r)V>j (r)^ (r) 



describes the dynamics of the atoms in the electronic 



state \j) at frequency 



Here ipj(r) and 



^t(r) 



arc 



the annihilation and creation operator, respectively for a 
bosonic atom at position r and in the electronic state \ j = 
1,2), obeying the commutation relation [^(r), ^(r')] — 
Sj,k^(r — r'). Parameter gj = 4irh 2 a s j /m is the strength 
of s-wave scattering between the atoms in state \j) , 
with a s j the corresponding s-wave scattering length [29jj . 
Hamiltonian term H\2 describes s-wave scattering be- 
tween an atom in electronic state |1) and an atom in 
electronic state |2), see [H, [HI, and will be negligible in 
the situations we consider. 

The Hamiltonian terms for the interactions between 
photons and atoms are given using normal ordering [35j j . 
They include a term describing coherent Raman coupling 
between the states |1) and |2), in which a photon in mode 
1 is scattered into mode 2 and vice versa. Raman tran- 
sitions follow from a pump-probe type of excitation and 
the corresponding Hamiltonian reads 



H, 



dr 



7(r)<4^!(r)^i (r)ai + H.c. 



(3) 



with the position-dependent coupling strength 7(1*) [3f 
For later convenience we denote by the quantity 



UJ12 — ^1 — ^2 



H — Hq + H em { + H int + iishift , (1) 

with Hq the Hamiltonian governing the atoms dynamics 
in absence of interaction with the lasers, H em f the Hamil- 
tonian for the free transverse electromagnetic (e.m.) 
field, while Hi nt describes the coupling between states |1) 
and 1 2) induced by coherent Raman scattering. Hamilto- 
nian ifshift includes the dynamical Stark shift due to off- 
resonant coupling between atoms and e.m. -field modes. 
These Hamiltonians are given using a second-quantized 
description of matter- and photon-fields. In detail 

H cm f = ^ hu x a\ax 

A 

gives the energy of the transverse electromagnetic field in 
free space (without the vacuum energy) where A labels 
a mode of the e.m. -field at wave vector k, polarization 
e _L k, and frequency Wk = c|k|, with c the velocity of 
light. The operators a\ and a\ annihilate and create, 
respectively, a photon in mode A, and obey the bosonic 
commutation relation [aA,a^,] = <5a,A'- 

The atomic Hamiltonian Hq is conveniently rewritten 



the difference between the frequencies of the two e.m.f.- 
modes. Moreover, we assume u>\2 ~ w a 2 — ^ai, i-e., the 
two lasers drive quasi-resonantly the Raman transition 
coupling the stable electronic states |1) and |2). The 
interaction with the pump and probe laser also induces a 
dynamical Stark shift, whose corresponding Hamiltonian 
reads (in normally-ordered form) 

ffshift = h J dr Tj(r)a^](r)^-(r)a j , (4) 
i=i 

where the parameter 7^ (r) has the dimension of an angu- 
lar frequency. It can be written as 7j(r) = 2|f2j(r)| 2 /A, 
with Oj(r) the Rabi frequency of the mode coupling 
the transition \j) — > |e) and A the detuning between 
the mode and the atomic transition frequencies (the 
Lamb shifts arc included in the frequencies of the atomic 
states). 

In the following we will assume that the pump and 
probe fields are traveling waves with wave vectors kj. In 
this case, the Raman- coupling strength 7(1") reads 

7 (r) = 7oe^ r (5) 
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with q = ki - k 2 and 70 = 2fti(r)fJ 2 (r)*/A. 

We remark that the outcoupling of atoms from two 
Bose-Einstein condensates in the setup of Refs. [3(| E[ 
was performed by Bragg scattering while we model the 
outcoupling process by Raman coupling to a different 
hyperfine state. This difference however does not affect 
the general results, as it will be shown in the following 



and where we have set u> a i = and g\ = g. In this repre- 
sentation the Hamiltonian describing the interaction with 
the lasers takes the form 



B. Spin-dependent potential 



We will now provide more details on the potential 
which confines the atoms. We will assume that the atoms 
in state |2) propagate freely, namely, potential V 2 (r) in 
Eq. @ is a constant which we set equal to zero. 

The atoms in state |1) are trapped by potential Vi(r). 
This potential confines the atoms in two spatially sepa- 
rated regions around the two potential minima, which are 
at distance d and, specifically, are localized at the points 
tl = —(d/2)x and tr — (d/2)x. The potential can be 
then decomposed in the sum V\ (r) = Vl (r) + Vr (r) , with 
Vj(r) the potential centered at rj and j = L, R. The two 
atomic clouds at each well are initially uncorrelated and 
there is no tunneling between the two spatial region. In 
principle, hence, the atoms in the two regions of space 
are distinguishable. 

It is useful to consider the partition tpi (r) = ^i(r) + 
iPr(t) where we denote by tp^r.t) — ipi(r,t)8(—z) and 
ijjji(r,t) = ipi(r,t)6(z) the field operators which do not 
vanish on the left and right region of space, respectively, 
such that [t/jj(r),ipl(r')] = S j:k S(r - r') for r,r' ^0 |3. 
Using these definitions we can write Hi = Hl + Hr with 



H 



{j—L,R} 



dr^(r) 



-n 2 v s 

2m 



+ V J (r)Up j (r) 



drV>t(r)V](r)^(r)^(r), (6) 



dr 



7(r)a 2 V4(r) (ip L (r) + Vfi(r)) 01 + H.c. 



(7) 

We note that our model is an extension of the one in 
Refs. [38|,|39| where the authors studied the outcoupling 
of atoms from a single BEC by means of classical Ra- 
man lasers. By considering the quantum dynamics of 
the interaction between photons and atoms we take into 
account the quantum fluctuations of the light field due 
to the scattering process, which are instead discarded in 
38, 39]. We will indeed show that these fluctuations give 
access to some correlation functions of the scattering sys- 
tem. 

C. The scattered field 



We now study the properties of the scattered photons 
by considering the photon field operators in Heisenberg 
picture. The Heisenberg equations of motion are deter- 
mined assuming that the coupling between matter and 
photons is sufficiently weak to be treated in second-order 
perturbation theory. In this limit the operators for the 
photonic modes read 



ai(i) = ai^e-^+^W) 



le 



-ituit 



dre" 



a 2 (t) = a 2 (0)e 



dr 7 (r) [< }t (r, r) + ^ (r, r)j ^ (r, r)a 2 (0) , (8) 
dr 7 (r)?/4 0)t (r, r) Uf (r, r) + (r, r)l a x (0) , (9) 



r 



where ipf\r,T) = exp(iH T/h)ipj (r, 0) exp(— \H q t/K). 
The physical origin of the individual terms on the Right- 
Hand Side (RHS) can be simply identified. The first 
term on the RHS of both equations is the free-field com- 
ponent. It is characterised by a time-dependent phase 
(t) , which is proportional to the atomic density in the 
electronic state \j). When the Born approximation can 
be performed, it is a density-dependent phase shift of the 



field mode, with the form 



4 



I 

0! (t) = J] / dT / dr Ti (r)Vf )f (r, r)^ 0) (r, r) , 

t 

0a(*) = |dr| dr 72 (r)v4 0)t (r,r)4 0) (r,r). 
o 

Given that all atoms are initially prepared in the elec- 
tronic state |1), term (j>\{t) gives rise to a frequency shift 
of field-mode 1. This shift is proportional to the density 
of the medium integrated over the path along which light 
propagates (3, Eoj. 

The second term on the RHS establishes a direct 
proportionality relation between the photonic and the 
matter-wave fields a 2 and ipi(r) (ai and ^(r)). It shows, 
in particular, that the source term of field a\ is the co- 
herent overlap of the field scattered from the right and 
from the left wells. In this shape the outcoupling pro- 
cess in this system resembles a beam-splitter operation. 
This property establishes an analogy with an interfero- 
metric setup, which we will exploit in order to perform 
homodyne detection of matter- wave field. Differing from 
a simple interferometer, however, atoms and photons are 
correlated by the scattering process. This important dif- 
ference is at the basis of the dynamics we observe. 

Before we discuss the signal at the photodetector, we 
will introduce a few approximations which will notably 
simplify the treatment. In first place we neglect atomic 
collisions between outcoupled atoms in state |2), since 
we assume that the gas of outcoupled atoms is at very 
low densities. This regime also permits us to neglect 
collisions between outcoupled atoms and trapped atoms 
in state |1) 0. 



III. PHOTON FLUX 

The quantity we analyse in this paper is the flux of 
photons of mode 2, namely, the rate of change of the 
photon number in the mode at frequency uj 2 . Formally, 
the photon flux is given by the equation 



F(t) 



d_ 

dT 



a\{t)a 2 (t)) 



(10) 



where the expectation value (•) is taken over the initial 
state of the atoms and the e.m. -field. The photon flux, 
integrated over the detection time, gives the integrated 
intensity of the field at the detector. It can be verified 
that 



dt 



F(t)K- / dr<v4(r,i)V> 2 M)) 



(11) 



where the proportionality factor is the mean number of 
photons in mode 1. This equality shows indeed that flux 
of scattered photons and of the corresponding outcoupled 
atoms carry the same information. 



For the specific setup we consider, the photon flux can 
be rewritten as the sum of two contributions, 



F = Fb + Fj 



(12) 



with 



F B (t) = TRc / dr / dr ' 

j=R,L 

1 

x J At' fir, t- r', t')G 00 (r, t; r', t 1 ) , (13a) 
o 

t 

Fj(t) = TRe J dr J dr' J dt'f(r, t; r', t') 

o 

x {G LR (r, t- r', t') + G RL {r, t; r', t')) [13b) 

where T = 2{a\a\)^Q is a scaling factor that is propor- 
tional to the number of photons in mode 1. Here, 

G jk (r,t;r',t') = (^t(r',t')Vfe(r,t)) , (14) 
f(r,t;r',t') = ^- r ')-^(t-t')] U 2 ( r ' ,t')4(r,t)" 



(15) 

where the initial state is the vacuum state of the e.m. 
field except for mode 1 and 2, which are assumed to be 
in coherent states with non-vanishing photon number, 
while all atoms are in internal state |1) and at equilibrium 
in the grand-canonical ensemble at temperature T. The 
atoms in state |2) propagate freely and do not undergo 
collisions since we assume that the density is very low. 
Therefore, Eq. (TT5j) can be cast in the form 



dk 



(2tt) 3 



f{r,t;v',t') = I —v 
with 



i[(q-k)-(r-r')-(wi2-Wk)(i-t')] 



(16) 



2m 



the recoil frequency and 



UJ 12 — Ld\2 — UJ a 2 

the two-photon detuning. 

In order to obtain explicit expressions for the corre- 
lation functions Gjk(r,t;r',t') it is convenient to work 
in the interaction picture with respect to the grand 
canonical ensemble Kq = Hq — X^=l J? Mj-Vj > where 
fij is the chemical potential of the atoms in either the 
right or left cloud and A/} = J dr ip^(r)ipj(r). The 
new atomic field operators are obtained from the ones 
in Heisenberg picture with respect to Hq by replacing 
ifjj(r,t) — > "0j( r ! £)e - "s wt . From now on we will denote 
by tpj (r, t) the atomic field operators in Heisenberg pic- 
ture with respect to Kq. 
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We now discuss how the considered setup can be used 
to measure the mean value of the atomic field operator 
by means of photons. We first note that the compo- 
nent Fb of the photon flux is the sum of the flux from 
each well, while the component Fi arises from the co- 
herent superposition of a photon (atom) scattered by 
the wells. We will denote Fb by "background contri- 
bution" and Fj by "interference contribution" . In ab- 
sence of initial correlations, this latter term is propor- 
tional to the product of the mean value of the field oper- 
ators (V'iCr, t)){tf)ii(T, t))* . Let us now consider the case 
in which, say, the left well confines weakly-interacting 
atoms forming a Bose-Einstein condensate. Under this 
assumption the atomic field operator can be written as 

M 

^(r,t) = e-*^*($i(r)+S^(r,t)) , (17) 

where $z,(r) is the macroscopic wave function which 
solves the Gross-Pitaevskii equation for the quantum gas 
in the left potential well (with chemical potential /xl). 
Then, when the order parameter of the left condensate is 
known, the interference term will deliver the mean value 
of the field operator in the right well. 

In the following we use Eq. (jTTJ) in the equations for the 
background contribution to the photon flux, Eq. (I13al) . 
and for the interference term, Eq. (| 13h>[) . In particular 
we will take 



$ i (r) = <^(r,t))=/ L (r)e i » 



(18) 



and assume that both /i(r) and ifL are real valued. 
Hence, /l(i") 2 is the density of condensed atoms, while 
the phase tp^ is assumed to be constant in space (hence 
discarding the possibility of superfluid currents in the 
Bose-Einstein condensate [Hj]). 



A. Background contribution 

The background contribution can be decomposed into 
the sum of the photon flux from the left and from the 
right well, Fb = Fl + Fr. An explicit form of Fl can 
be provided since we make a specific assumption on the 
quantum state of the gas trapped in the left well. For this 
purpose, we consider the integrand Gll{t, t; r', t') in Fl 
and observe that it can be splitted into two terms: 



G LL {v,t-v\t' 



= c WL(t-t')/h 



[f L (r) 2 + 6G LL (r,t;r',t')] , 



accounts for the contribution of the noncondensed atoms 
to the photon flux. Here, 

oo 



(19) 



where 



5G LL {v,t-v\t') = (5^ L (r,t)5^ L r',t')) 



2^ 



-iw(t-t') 



W^^tM,^) (20) 



n (21) 
is the spectral density |42j, with N Q (u) = [e /3ri " - l]" 1 . 

Similar expressions can be found for the background con- 
tribution of the quantum gas in the right well. 

Insight on these equations can be gained by assuming 
that the gas is homogeneous. This assumption allows 
us to write the spectral density as A s ^ s ^ (r, r', cj) = 

■A-g^.g^X (l r—r 'li We then inspect the equations for the 
flux using the momentum representation for the atomic 
field operator. The terms giving the contribution to the 
photon flux from the left cloud can be cast in the form 



F L (t) = 27rrRe 
SF L (t) 



dk 

(2tt) 3 



2tt 



l/iOOI^n-Wk+q), (22a) 
dk 



(2tt) 3 

X N Q {u))A HLH i (k, -w)<5*(w + n - w k+q ) 



(22b) 



where A 



8il> R 5il)' R 



(k, uj) is the spectral weight function, 



which is the Fourier transform of A 



SlpL Sip 



(r,ui) and gives 
the strength of the collective excitations at wavevector 
k and frequency w, and is weighted by the Bose func- 
tion Nq(u)) [42]. The other parameters are the vol- 
ume V in which the atoms of the left well are confined, 
the Fourier transform of the macroscopic wave function 
/i(k) = / dre _lk ' r /i(r), and the frequency 



n 



"12 



h 



(23) 



giving the effective detuning of the laser from the collec- 
tive transition. Moreover, in Eqs. (|22l) we have intro- 
duced the quantity 



(24) 



that is proportional to the diffraction function 5 (u>) — 
sm(u>t/2)/(iruj), enforcing energy conservation for long 
interaction times t. In particular, 5 t {uA — > 8(us) for 
t — > oo, with 5(uj) Dirac-delta function [44|. Definition 
([Ml) contains a time-dependent phase and highlights that 
the corresponding factor tends to unity in the limit in 
which energy conservation applies. We will keep this 
time-dependent phase in the equations for the photon 
flux since we will also consider intermediate times. 

The integrals in Eqs. (|2"2"|) runs over all values of the 
atomic momentum k. The integrands are the product 
of the momentum distribution at k' = q — k and of 
the diffraction function: The first accounts for the 
effect of the photon recoil hq on the atomic distribution 
due to the outcoupling process while the diffraction 
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function imposes energy conservation in the scattering 
process. Equations (f2"2"b) and (|22b ) thus show that the 
contributions to the flux come from the total number 
of condensed and noncondensed atoms, respectively, 
which fulfill energy and momentum conservation of the 
scattering process. 

We remark that according to Eq. (ITT|) the flux evalu- 
ated in Eqs. (|22|) gives also the corresponding component 
of the atomic flux. The latter agrees with the correspond- 
ing expressions derived in Ref. [HJ for the atomic flux 
which is outcoupled from a single BEC by classical fields. 

B. Interference contribution 

In order to determine the interference contribution to 
the photon flux, Fj(t) in Eq. (|13b[) . one needs the ex- 
plicit form of the correlation functions Gl R {t, t; r', t'), 
Crii£(r,t;r',i'). For the considered setup, however, one 
can already make general statements. In absence of ini- 
tial correlations, in fact, they are the product of the mean 
value of the field operators in each well, and thus take the 
form 

G LR {v,t;v' ,t') = ^ Lt -** e)/A (1>l(T,t)){<l>RpJ)) 

(25) 



while G RL (r,t;r',t') = [G LR (r', f; r, t)]* . 



Equation (|25|) shows, as it is expected, that these cor- 
relation functions are proportional to the mean value of 
atomic field operator (ip R (r,t)). When this is zero the 
interference contribution vanishes. Otherwise, the am- 
plitude of this term is proportional to the condensate 
fraction of the quantum gas in the right well. We re- 
mark that the noncondensed atoms in the left well do 
not contribute to the signal because (i) there are no ini- 
tial correlations between the atoms in the left and right 
wells and (ii) the mean value (StpL(^,t)) = 0. 



In the shape of Eq. (|2"5"|) the analogy with homodyne 
detection, as it is performed with light fields, can be 
drawn. The condensate in the left well here plays the 
role of the local oscillator (28|. Identifying the analog 
of the phase of the local oscillator is, however, a more 
delicate issue that deserves some more analysis. For this 
purpose we first assume that the mean value (ip R (r,t)) 
can be evaluated within a mean-field approach, such that 
(ip R (r,t)) — f R (r)e l(fR , with f R (r) real-valued function 
and <p R real constant. The integral in Eq. (113b[) can then 
be cast in the form 



Fi(t) = 27rrRe 



dk 

J2nf 



3 i(6jLit— <plr) 



/ L (k)*/ fl (k)^(ft - c k+q ) + e- i ^-^")/ fl (k)*/ L (k)5*(fi -5n- c k+q )l , (26) 



with 



Discussion 



8fi = (n L - (x R )/h, (27) 

and 



fL R = <Pl- <Pr ■ (28) 

Equation (|26p indicates that the interference contribu- 
tion is an oscillating signal. The phase of the oscillation, 
however, depends on time and oscillates with a frequency 
determined by the difference 5fi between the chemical po- 
tentials while the phase offset is determined by the rel- 
ative phase between the two (quasi-)condensates. This 
relative phase will be defined only for a single experi- 
mental run, being the two quantum gases independent 



1. Which-way information and quantum erasers 

Some remarks on the oscillating behaviour of the in- 
terference term are now in order. For the example con- 
sidered in Eq. (l2l)l) the photon flux oscillates in time 
with frequency 5[i. This oscillation is observed when 
one chooses a time-window at the detector At such that 
S/iAt <C 1. The interference arises from the overlap of the 
atomic beams from both condensates. It is analogous to 
the interference between two lasers at different frequen- 
cies: interference fringes are observed provided that the 
time window of the detector is sufficiently small so not 
to resolve the detuning between the lasers [HI [52| . As in 
the case of two independent lasers, there is no well de- 
fined phase offset: the relative phase is defined only for 
a single experimental run while the average over a sta- 
tistically significant number of runs gives no interference 
pattern 0, El H, H . 
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Differing from the situation of two laser beams [51], 
however, photon scattering here creates correlations 
among the scattered photon, the corresponding outcou- 
pled atom, and the quantum gas. This correlation is a 
form of a which-way information, that in general washes 
out any interference in the photon flux j53| and can be 
considered as a form of photon-atoms entanglement [54j . 
Oscillations in the photon flux, and hence interference, 
can be recovered if the following two conditions are ful- 
filled: (i) the beam of atoms outcoupled from one well 
overlap spatially with the wave function of the atoms in 
the second trap and (ii) if there are either initial correla- 
tions between the atoms in the two wells or if they both 
possess a non-vanishing condensate fraction. 

In order to better understand these conditions we 
first notice that the outcoupled atoms are distinguish- 
able from the trapped atoms since they are in different 
quantum states. Nevertheless, the Raman beams can, in 
principle, outcouple one atom from one trap and then 
load the same atom into the second trap. This effective 
coupling is at the basis of the analogy with a Josephson 
Junction that was drawn in Refs. |30L l3l| . and is the es- 
sential element leading to interference. In these terms, 
it corresponds to a quantum eraser [32j : when it is ful- 
filled, in fact, the outcoupled atom becomes disentan- 
gled with the scattered photon and with the scattering 
system. Condition (ii) consists in assuming that there 
exist classical spatial correlations between the wells. It 
is important since we assume that at t = there are no 
quantum correlations between the two scattering systems 



we introduced the recoil velocity 

ftq 



(31) 



2. Onset of the time- dependent oscillations 

We now analyse the onset of oscillations in time. We 
have argued that the photon flux starts oscillating af- 
ter a finite time has elapsed from the beginning of the 
experiment. [3(| Hl|. This can be also seen in our the- 
ory by performing the integral over k in Eq. (|26l) . In 
App. we show that Fj can be recast in the form 
Fi{t) = F L ^ R {t) + F R ^ L {t), with 

t 

Fl-m® ps rRee 1 ^*-^) J dre i( ^-" )r /i LJt (d, r) (29a) 

o 

t 

F R ^ L {t) « FRee-' 1 ^-^ J dTe i{ ^- n+5f * )T h RL (-d, r) . 

o 

(29b) 
Here. 



that is acquired by the outcoupled atom by scattering 
the photon. This overlap vanishes at time r = 0, i.e., 
when the outcoupling lasers are switched on (recall that 
the two clouds initially do not overlap). We note that 
the overlap integral is zero at all times if d and q are 
orthogonal. It may be finite and reach a maximum after 
a certain time when d and q are parallel, say, pointing 
along the positive x axis. In this case the component 
F^ R (t) may not vanish and can be interpreted as the 
contribution to the interference flux Fj(t) from the out- 
coupled atoms that propagate from left to right (for the 
term F R ^L{t) this is just opposite). Interference will 
be observed provided that a sufficiently long time has 
elapsed to warrant overlap. This corresponds to times 
t > t r with 



tr. 



where d c g — d — (£l + £r)/2 is the effective distance of 
the two systems taking into account their width £j. In 
other words for times t > t c the which-way information 
has been erased and oscillations in the photon flux can 
be observed. 

In interferometric setups the amount of visibility and 
which-way information are related by an inequality [53|. 
It is therefore useful to determine a visibility of the os- 
cillating signal because it contains information on the 
properties of the scattering systems. A visibility can be 
defined for sufficiently long times by averaging over sev- 
eral oscillating periods after the instant t c . In this case it 
will be proportional to the amplitude of the oscillations 
of the photon flux, hence, to the product of the conden- 
sate fractions of both systems. We remark, once again, 
that oscillations can be observed only in a single exper- 
imental run while they will disappear after performing 
an ensemble average. Therefore, this behaviour can only 
be measured in systems whose properties are not deeply 
modified by the outcoupling of atoms. When this is not 
verified, the presence of a condensed fraction in the right 
system can be revealed by performing an ensemble av- 
erage over a sufficient large number of experiments in 
which the signal is monitored for sufficiently short times, 
warranting that the properties of the quantum gas have 
not been significantly modified, and taking the statistical 
distribution of the intensity of the photon flux at a given 
instant of time. An amplitude can be extracted from 
the width of this distribution by taking into account the 
finite width of the diffraction function. 



hji(d, t)= dr fj (r + d - u q r) / z (r) 



(30) 



is the overlap between the left and the right condensate, 
with one being displaced by the amount u q r — d. Here, 



IV. HOMODYNE DETECTION OF A 
QUANTUM GAS 



The theory presented so far will be now applied to 
some specific examples. The idea is to use the setup in 
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Fig. [I] to determine the mean value of the field operator 
of a quantum gas by using a Bose-Einstein condensate at 
known temperature as reference system. Such setup is 
a matter- wave analogon of homodyne detection in quan- 
tum optics. The individual elements can be so identified: 
the BEC acts as a local oscillator, the outcoupling pro- 
cedure as beam splitter, the relative phase can be varied 
by changing the interwell distance. The information on 
the atomic gas is carried by both scattered photons and 
outcoupled atoms: homodyne detection of the scattered 
field, hence, allows one to determine the mean value of 
the quadrature of the atomic field operator. 

A. Interference between two Bose-Einstein 
Condensates at T = 

We first discuss the case in which two BEC are trapped 
in the left and right well, respectively, and are both illu- 
minated by the pump and probe beams. The outcoupled 
atoms from both condensates propagate along the direc- 
tion determined by the vector joining the minima of both 
wells. The two atomic beams spatially overlap after the 
time t c has elapsed. The scattered photons are revealed 
at a detector in the far-field. This setup has been realized 
in the experiment presented in Refs. [30l.l3l| where time- 
dependent oscillations in the atom and photon flux were 
measured. We remark that the outcoupling process in the 
setup of Refs. [3(1 Hl| was performed by Bragg scattering 
while we model it by Raman scattering into a different 
hyperfine state. This difference does not affect the re- 
sults: While in Refs. [13, [HJ the atoms are transferred 
into a momentum state that freely propagates, here the 
atoms are transferred into a hyperfine state that is not 
trapped. In both cases, collisions between the outcoupled 
beam and the trapped atoms can be neglected. 

We assume two BEC with equal number of atoms Nc 
and temperature T, that are confined either in the left 
or right well (The assumption of equal number of atoms 



is a convenient theoretical choice, but no substantial re- 
striction to the results presented in the following). The 
wells are described by the potential 

fy=L,H}(r) = V(r + Sj,lAV (32) 

with V(r) = im(kiji 2 + uj 2 y 2 + uj 2 z 2 ), and AV denotes 
the constant offset between the two traps. For simplicity 
we take that both BEC are at zero temperature, T = 0, 
and the atoms weakly interact, such that the contribution 
of the noncondensed atoms to the photon flux is small 
and can be neglected. In this limit the chemical potential 
of both condensates is equal and given by /i(0). With 
definition (j3"2")l then /«r = //(0) and fj,£ = /i(0) + AV. 

The photon flux is evaluated using the Thomas-Fermi 
approximation for the condensate wave functions [29j , 

/ J -(r) = [( M (0)-y(r-^))/ g ]V2 i (33 ) 

where /x(0) = (l5Nca s /a) 2 ^ 5 huj/2 is the chemical po- 
tential at zero temperature and a = y/h/(mui) is the 
size of the ground state of a harmonic oscillator with fre- 
quency Q — (uj x ujyUj z ) 1 ^ 3 . The condensate macroscopic 

wave function has size rf^_ x y z y = ^2^(0) /( m <-° 2 )- Its 
Fourier transform reads /z,(k) = e lk ' d / 2 /o(k) (for the 
right condensate /n(k) = e _lkd / 2 / (k)) where /o(k) is 
the Fourier transform of the macroscopic wave function 
centered at the origin and is real valued. In particular 

/o(k) = k 5 — ' I 34 ) 



where «o = y 15n 3 Ncr x r v r\ /2 is a scalar and J^iPo) 
the Bessel function of second order for the variable po 
defined as p 2 = k 2 x r { x )2 + k 2 y r ( y )2 + k 2 7 rf )2 , see [29]. The 
integral for the interference contribution to the photon 
flux in Eq. (|26[) can be now cast in the form 



Fi(t) = 2ttkIt 



dk |J 2 (po)| 5 
(2tt) 3 



Re 



e i(«M*-V-d-*»£Jt)Jf( n _ Wk+q ) + e -i(S^-k-d- VLR )^t^ n _ Sfl _ Wk+q ) 



(35) 



This equation can be simplified taking that both d and q point along the positive ir-direction. Following the derivation 
in App. |B]we find that the total flux can be approximated by the expression 



F(t) « 2-KTNcKr, 



1 + Vq COS ( Sflt - lf LR + ((jjq - CI) — ) ( t 



d-2r,, 



(36) 



r 



where 

K m = (K(u}q - O) + X(a; q - ft + Sfi))/2 
and K(x) — y/TJ (2na 2 ) exp[— x 2 / (2a 2 )] is a Gaussian of 



width a 2 = 2.5(vq/r x ) 2 . Equation (|36|) shows that the 
photon flux starts oscillating for times t > t c , with t c — 
(d — 2r x )/vq, namely, when the outcoupled atoms from 
one condensate have reached the second one. The time- 
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dependent oscillations have frequency 5fi and amplitude 

= M f371 

if ( Wq - 0) + Kiuiq -Sl + 6fi)' V ; 

The expression in Eq. (|5f?)) has been obtained neglecting 
the momentum dependence of the trapped clouds (this 
approximation is verified in App.[B]) and the contribution 
of the noncondensed fraction. We note that, although 
oscillations as a function of time are observed provided 
that ^ 0, nevertheless their visibility (which corre- 
sponds to Vo) is always smaller than unity, Vo < 1 for 
5fi =/= 0. This can be understood considering that a scat- 
tered photon carries information about from which cloud 
it was scattered because of energy conservation in the 
scattering process. 

We now turn our attention to the phase offset charac- 
terising the oscillations of the photon flux as a function 
of time in Eq. (f5B|) . This phase offset is here given by 
the quantity O — <plr, where <plr is the relative phase 
between the two BEC, which can take any value (we refer 
to the discussion in Sec. IIIICp . and by the quantity 

6 = — (w q - Q + 8/j,) , (38) 
v q 

which is the phase the outcoupled atoms accumulate 
when travelling from left to the right well (The corre- 
sponding phase offset for the outcoupled atoms travelling 
from the right to the left well is given in App. [B]). This 
phase can be varied by either tuning the laser frequency, 
and thus f£, or changing the distance between the min- 
ima of the two wells. The phase offset in Eq. (|38[) agrees 
with the one obtained in Ref. [3l[ that was derived from 
a phenomenological model [56j . 

We also note that our theory correctly predicts the 
oscillation period observed in the experiment and also 
reproduces the behaviour that the interference current 
needs some time to build up. The linear response treat- 
ment we apply, however, does not predict any decrease in 
the visibility of the interference pattern with time. This 
in marked contrast with the experimental results. Possi- 
ble reasons for the decrease in visibility in the experimen- 
tal data are depletion of the condensate and heating of 
the sample, which could be due to spontaneous Rayleigh 
scattering events [3(| ■ These effects are not taken into 
account in our model, but could be introduced by means 
of quantum Langevin equations using a formalism similar 
to the one developed in Ref. [35| . 

We finally discuss the dependence of the photon flux 
on time and on the angle of emission of the outcoupled 
atoms, thus for geometries where the direction of emis- 
sion q of the outcoupled atoms forms an angle a with 
the vector d. For simplicity we assume that both vectors 
lie in the x — y plane. The calculation is performed by 
numerically integrating Eq. ([35]) . 

Figure [2] displays the photon flux F(t) as a function 
of time and of the emission angle a. Oscillations as 
a function of time are observed for values of the angle 




FIG. 2: (color online) Photon flux F(t) (in units of the back- 
ground contribution Fs(t)) as a function of time (in units 
of 8fj,/2iv) and of the angle a at which the atoms are emit- 
ted. The photon flux is computed by numerically integrating 
Eq. (|35l) . The condensates are composed by N = 10 6 sodium 
atoms in a spherical harmomic trap of frequency w = 325Hz. 
The scattering length is a s — 55ao with ao Bohr radius. The 
other parameters are « q = 6^, <5/io = 27rl0 3 Hz, d = 5r x , 
£1 = cj q , q = and i^lr = 0. 

about a = 0, 7r corresponding to the atoms propagat- 
ing in the direction parallel to ±d. The oscillation pe- 
riod is 2tt/5[a and is independent of the angle of emis- 
sion. The oscillations disappear for angles in the interval 
11 /& ^$ ol < 7 ir/ 8: the photon flux is here solely given by 
the background contribution Fs(t). For these angles, in 
fact, at all times there is no spatial overlap between the 
wavefunction of the outcoupled atoms and the wavefunc- 
tion of the trapped atoms in the other well. Indeed, from 
a simple geometric argument one finds that the overlap 
vanishes for angle a > arctan w ir/8 and a < tt — tt/8. 
This implies that the setup must be so constructed that 
the atoms outcoupled from one well could in principle be 
transferred, by a Raman process, into the second well. 
This property is the key element on which the analogy 
to a Josephson Junction has been drawn |3lj . It is also 
basically the way in which a quantum eraser is realised in 
this setup. We refer the reader to the discussion on the 
properties of this quantum interference process in Sec. 

mrci 



B. Interference between two Bose-Einstein 
condensates at different temperatures 

We now show how the setup of Ref. 30] could be used 
to determine the condensate fraction of a Bose-Einstein 
condensate at a different temperature, using as reference 
a second BEC at known temperature. For this purpose, 
we make the same assumptions as in the previous section, 
with the difference that while the left gas is a BEC at 
temperature T = 0, the gas in the right well is a BEC 
at temperature T. We further assume that the Thomas- 
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Fermi approximation can be performed also for the right 
condensate. Hence, the chemical potential of the second 
condensate can be written as fj,(T) = ^(0)(1-T 3 /T 3 ) 2 / 5 , 
with T c the critical temperature for the noninteracting 
gas in a harmonic trap, while the size of the macroscopic 



wave function of the right condensate scales with r\ — 

rf\l - T 3 /^ 3 ) 1 / 5 , see 0. Using these relations the 
integral for the interference contribution to the photon 
flux can be cast in the form 



Fj(t) = 2-KTnlnciT) J ^ 



dk | J 2 (po)^2(pt)| 



xRe 



(2tt) 3 



Pt 



) + e -i(« M (Dt-k.d-#, ia )^(Q _ 5^(r) - Wk+q ) 



(39) 



wither = (1— T 3 /T^Y^pq. The integral depends on the 
temperature both through a scaling factor as well as the 
function Ji{j>t)Ivt m the integrand, while the phases 
depend on the temperature of the second BEC via the 
chemical potential of the right BEC which enters in the 
quantity 



Sfi(T) = S^(0) -F{N,T)/h, 



(40) 



with F{N, T) = /i(T) - fi(0). We note that £>(t) oscil- 
lates as a function of time with both frequency and am- 
plitude which depend on the temperature of the second 
condensate. For T <C T Cl in particular, the integral in 
Eq. d39]) delivers the factor l/yJn c (T), and thus F/(f) cx 
\/nc(T), where n c {T) = N C (T)/N = 1 - {T/T c f is the 
condensate fraction in the right trap. 

Figure (j3£i) displays Fj (t) as function of time and tem- 
perature T of the right condensate when the trap is a 
spherical harmonic oscillator. The oscillations are visible 
for times t > t c . The oscillation frequency depends on 



C(T) = [max(F 7 (i)) - min(F/(t))] 



(41) 



is displayed as a function of T in units of C(0). The 
red dashed curve represents the squared root of the con- 
densate fraction, y/nc(T), and is shown for comparison. 
Function C(T) decreases monotonically as the tempera- 
ture increases from T = 0, and vanishes at the critical 
temperature T = (where the superscript i indicates 
that we take into account the effect of the interactions 
and T}j is determined from the condition nc(T*) = 0) 
581. 



stance Refs. 0, HH, |6(J. Measurements of the first- 
order correlation function have been performed on Bosc- 
Einstein Condensates in Refs. [6ll |62||. 

In the following we will show how an extension of our 
previously considered setup may allow one to determine 
the spatial first-order correlation function of a quantum 
gas by the scattered photons. In this case, the lasers shall 
illuminate two spatially separated regions of a quantum 
gas confined in a single well potential. More specifically, 
the spatial dependence of the laser-atom interaction in 
Eq. [3] will be characterized by the Gaussian envelope 

| 7 ,(r)|= 70 exp(-(r-r j ) 2 /Ar 2 )/(^Ar) 1 / 2 , 

with width Ar <C |d| between the two regions (We re- 
mark that the excitation could be realized with subwave- 
lcngth resolution d, 0, [fil])- The corresponding setup is 
shown in Fig. [4] For this setup the photon flux reads 



T, as one can clearly observe from the figure 57]. The 
dependence of the amplitude of the oscillation on the 
temperature is visible in Fig. and is singled out in 
Fig. ((3Jd) where the amplitude evaluated at times t > t c , where 



F(t) 



T 2^ Re / dt' I dr / dr' 

j,k=L,R 







x 7j (r)7fc (r')/(r, t; r', i')G (1) (r, t; r', t>) (42) 



G«(r ! t;r',t') = <^(r,^i(r / ,t / )) 



(43) 



is the first-order correlation function and f(r,t;r',t') 
is defined in Eq. (|15|) . Let the pump-probe excitation 
be a pulse of mean duration t such that Qt ^> 1, but 
uj a t <C 1, with Lo a the typical frequency characterising 
the excitation spectrum. In this limit we can approxi- 
mate G^\r,t;r',t') ~ GW(r,0;r',0) in Eq. (g2J . For 
convenience we denote by G^\r; r') = G^ 1 ' (r, 0; r', 0) the 
spatial correlation function. For Ar sufficiently small, so 
that it can be approximated with a (5-likc excitation, the 
photon flux can be recast in the form 



V. FIRST-ORDER CORRELATION FUNCTION F(t) ~ K (l + Re |e iq d G (1) (d/2; -d/2)| /n (d/2] 



Several methods that are based on time-of-flight tech- 
niques for determining the first-order correlation func- 
tions have been proposed in the literature, see for in- 



(44) 

where if is a constant which is determined by the details 
of the excitation scheme and no(d/2) = G^- 1 ' (d/2; d/2) 





T 

(b) 

FIG. 3: (color online) (a) Interference term of the Photon flux, 
Fi(i), Eq. (|39[) . (in units of the background current obtained 
when both condensates are at zero temperature Fb{T — 0)), 
as a function of time (in units of 2ir/8fj,) and temperature T 
(in units of critical temperature T^). (b) Amplitude C(T), 
Eq. (|41|) (in units of C(0)) as a function of the temperature 
T (in units of The red dashed line corresponds to the 
squared root of the condensate fraction in the right BEC, 
y/nc{T). The other parameters are as in Fig. [2] 



is the density at r = ±d/2 (assuming the system has 
reflection symmetry about r = 0). Therefore, the photon 
flux exhibits oscillations with a visibility which is deter- 
mined by the spatial first-order correlation function. By 
varying the scattering wave vector q one would thus mea- 
sure the first-order correlation function as a function of d. 
Realistic excitation schemes are of course characterized 
by finite spatial resolution Ar. This results in averaging 
the correlation function in Eq. (|44|) over the finite size of 
the illuminated region and hence to a diminution of the 
contrast [55j |. 

This setup could be extended to measure time- 
dependent correlation function by applying a pair of laser 
pulses: Assuming that the photon scattered after the first 
pulse can interfere with the photon scattered after the 



FIG. 4: Setup for measuring the first order correlation func- 
tion of an ultracold atomic gas. The pump probe lasers are 
tailored in such a way that atoms are only outcoupled from 
the two shaded regions around ri ; 2- Measurement of the pho- 
ton flux of one beam as a function of q allows one to determine 
the first order correlation function G^(ri;r2) of the atomic 
gas. 

second pulse, then the photon flux at the detector will ex- 
hibit oscillations whose amplitude is proportional to the 
first-order correlation function Eq.([43|). A possible real- 
ization could use a mirror placed in front of the quantum 
system as realized in [64l |. 



VI. CONCLUSIONS 

In this work we have discussed a setup which allows 
one to measure the condensate fraction and the first- 
order correlation function of a quantum gas by means 
of photo-detection. The photons are scattered by the 
quantum gas in a pump-probe type of excitation, such 
that the scattered photon is associated with an outcou- 
pled atom with which it is entangled. In addition, pho- 
ton and atom are correlated with the quantum gas. This 
correlation is detected in the photon flux, provided that 
certain conditions are fulfilled which we have identified 
and discussed. 

Our analysis is based on the impulse approximation 
[45j corresponding to neglecting the back action of the 
scattering process on the quantum gas. It is therefore 
valid for short time transients. From the point of view 
of collecting a sufficient statistics, hence, time-of-flight 
techniques are a more convenient tool than in-situ mea- 
surements by photodetection. Nevertheless, one could 
consider to modify existing techniques, such as the one 
successfully demonstrated in Refs. [g, 0], to access to the 
same kind of information that the setup here discussed 
provides. 
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An interesting outlook is to identify a setup, along 
the lines of the proposal Ref. [65j . which permits one 
to perform a quantum-non-demolition measurement of 
any correlations function of the external degrees of free- 
dom of atomic gases, and more in general, which realizes 
quantum-state transfer between matter and light. This 
would open several interesting perspectives for quantum 
communications [66j . 



We now show that the term F^nit) (F R ^L{t)) is the 
contribution due to the atoms which are outcoupled and 
propagate from the left to the right (right to the left). 
For this purpose we perform the k integral in Eq. (|Aip 
and obtain 
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Appendix A 



We recast the interference term, Eq. (|13b|) . in the form 
Fi = F L ^ R {t) + F R ^ L (t), with 



-i(n-w k )(t-t') 



fU*(t) = TRee^Jdt'J ^ 
o 

x J drdr'e i ^ k ^ r - r ')$ L (r')*$fl(r),(Ala) 

t 



Ldk — 5fi)(t— t') 



J dr dr'e i(q - k) ' {r - r ' ] t> R (r')* $ L (r) 



(Alb) 



F L ^ R (t) = rReeW"*-***) / dr 



2irhT 



dr / dr'e i(q - (r - r,) - OT) exp 



i m(r' — r) s 
h 27 



h(r>)f R (r), (A2) 



where <plr is the relative phase of the macroscopic wavefunctions defined in Eq. (|2"8)> . With the change of variables 
f = r r', R = (r + r')/2, we can rewrite Eq. (|A2[) as 



F L ^ R (t) = TRee^-^ J dr (j^A ' e i( "<*-^ T / <lr / dr 



exp 



i m I _ hq 
H 2r \ m 



/ i (R-£)/ fl (R+|)(A3) 



r 



The exponential in the integral over r oscillates very the main contribution to the integral over r comes from 
fast with respect to the wave functions fj(r). Therefore 
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fo = — t, where the term in the exponential vanishes. 
By means of the saddle-point approximation we take the 
wave functions at the point r out of the integral. Per- 
forming the f integration using the Fresnel integral (28j 



oo 

J "VH 



(A4) 



we obtain 



o 

x j drf L (r + d - « q r) /n(r) , (A5) 

where v q is the recoil velocity, Eq. ([3~Tj) . The calculation 
that leads to Eq. (|A5I) is exact in the limit of homo- 
geneous atomic systems. In such case the momentum 
distribution of the condensate fraction is a Dirac-delta 
function at zero momentum and all outcoupled atoms 
have exactly the same momentum hq. The condensates 
we consider are confined by an external potential and 
have a finite extension in space which leads to a certain 
width dp in their momentum distribution, and thus to a 
spread in the momentum of the outcoupled atoms around 
the mean value hq. By taking the value of the atomic 
wave functions at the point of the stationary phase of 
the exponential in Eq. (|A3I) one neglects this momentum 
width: The saddle-point approximation is thus applica- 
ble if h\q\ » Sp. Equation (|A5[) agrees with the corre- 
sponding expression in Eq. (|29a ). For completeness we 
also give the contributions of the macroscopic wavefunc- 
tions to the background contribution, Eq. (113a[) . making 
the same approximations as in Eqs. (|A5|) : 



z 

F L (t) » TRe J dre 1 ^-^ J dr/ L (r)/ L (r - u q r) , 
o 

t 

F R (t) « TRe J dTe iK-n+^)r J dr/ fl (r)/ fl (r- « q r) . 
o 

(A6) 

Appendix B 



We derive here approximate expressions of the Raman 
scattering rate for the experimental setup of [30| . which 
lead to Eq. (131)1) . Using Eq. (j3"3")l for the condensate wave- 
functions in Eq. (|A6[) we find 



F L {t) 



2irTn(0)r x R 3 

9lVq 



*<(0) 



dz cos 



(a; q - n)- 



G(z), 
(Bl) 



with 



z < (x) = min(2, — — ) . 



(B2) 



The length R = (r^r,^) 1 / 3 determines the typical size 
of the condensates and can be written as 67a 



= i,Va f^l 



R = 15 



1/5 



(B3) 



The function G(z) in Eq. (|Bip reads 



4/15 




FIG. 5: (color online) Comparison of the integral G(z) as 
given in Eq. (|B4|) (black solid line) with the Gaussian fit 
Eq. (|B5[) (red dashed line). 



G(z) = 2 I dx I rdr(l-x 2 -r 2 ^ 1/2 



z/2 



(1 - (x - z) 2 - r 2 ) 



2\l/2 



(B4) 



A numerical evaluation of G(z) is shown in Fig. ([3]), and 
is here compared with a Gaussian of the form 



G( Z ) « -e 



-1.252T 



(B5) 



showing that this function provides a good approxima- 
tion of Eq. (jB4|) . Using Eq. (|B"5]l in Eq. (|B"Tj) one gets 

z<(0) 



Fi(t) IW c i y dzcos[(w q -ft)zi ]e~ 
n 



(B6) 



with to = rx/vq. For times i > t c , with f c ~ (d—2r x )/vq, 
such that 2<(0) = z < (d) = 2 we find from Eq. (IB6|) 



A 

F L {t) « rt A c y dzcos 



(Wq - ^) — 



-1.25a;'' 



V 57T 



(B7) 
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When performing the integral in Eq. (|B7[) we neglected 
the imaginary part in the Error-functions which one gets 
from the exact integration. This is a good approxima- 
tion for |(cj q — f2)io| < 5, as can be checked by numerical 
evaluation. For values |(aj q — f2)io| > 5 the exponential 
in Eq. (|B7j) is negligible compared to the resonant case 
w q = Q: The oscillating tails of the Error- functions given 
by their imaginary parts will only play a role for parame- 
ters where the outcoupling efficiency vanishes and which 
are thus not relevant to our treatment. Using the same 
argumentation for the other contributions to the photon 
flux we find 



with 



K(x) = \hr-e-Tr- 

V 07T 

F e (q,t) = e(q)&(t 



d - 2r a 



(B9) 
(BIO) 



From Eqs. (IB8|) one then obtains Eq. (j3"fj)) for q > 0. 
The phase difference between FL^.ji(t) and F R ^L(t) as 
measured in [31| is obtained from Eqs. (|B8b .d) and reads 



F L (t) 
F R (t) 



■KTN C K(u}q - Q) , 

ttTNcK^ -Cl + Sfi), 
2nTN c K(u; (l -fl)Fe(q,t) 



X cos 



s fit - ipLR + (w q — fi) — 



(B8a) 
(B8b) 



c) 



F R ^ L {t) » 2ttTN c K(lu ci - n + S f i)F e (-q,t) 



6 (it - tfi LR - (ujq — Q + SfJ,)- 



(B8d) 



6= — (2fi - 2w q - <5/i) . 



(BH) 



Figure © displays the total photon flux computed 
from Eq. (f36|) (dashed line) and the one calculated from 
Eq. (|35|) (solid line) for the experimental parameters of 
30]. Quantitative agreement between the two solutions 
is found, showing that neglecting the initial momentum 
distribution is well justified for the considered parame- 
ters. 
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